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Abstract 

We present a light-front determination of the pionic contribution to the nucleon self-energy, Tl n , 
to second-order in pion-baryon coupling constants that allows the pion-nucleon vertex function to 
be treated in a model- independent manner constrained by experiment. The pion mass u dependence 
of St,- is consistent with chiral perturbation theory results for small values of fi and is also linearly 
dependent on \x for larger values, in accord with the results of lattice QCD calculations. The derivative 
of Ti n with respect to fi 2 yields the dominant contribution to the pion content, which is consistent 
with the d — u difference observed experimentally in the violation of the Gottfried sum rule. 
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Understanding the pion and its interaction with and amongst nucleons is a necessary step 
in learning how QCD describes the interaction and existence of atomic nuclei. As a nearly 
massless excitation of the QCD vacuum with pseudoscalar quantum numbers, the pion plays a 
central role in particle and nuclear physics as a harbinger of spontaneous symmetry breaking. 
The pion is associated with large distance structure of the nucleon [JJ [2] and the longest 
ranged component of the nucleon-nucleon force [3]. In lattice QCD calculations the nucleon 
mass depends on an input value of the quark mass, which generates a pion mass //, and 
extrapolation formulae depending on /i are typically used [IH7] (see the review [5].) In addition, 
the pion cloud plays an important role in deep inelastic scattering on the nucleon, especially 
in understanding the violation of the Gottfried sum rule [HI [TO] . 

Phenomenological calculations of pion-nucleon interactions are beset with uncertainties re- 
lated to the dependence of the vertex function on momentum transfer and on the possible 
dependence upon the virtuality (difference between the square of the four-momentum and 
mass squared) of any intermediate nucleon or baryon. Moreover, modern treatments of spin 
3/2 baryons such as the A (baryon excitation of lowest mass) within the Rarita-Schwinger 
(RS) |TT] formalism have been problematic as discussed in [12]. The pathologies of the nNA 
coupling have long been known [T3T[T7] . The aim of the present letter is to develop and apply 
a method that is free of those ambiguities. 

As a specific example, consider the role of the pion cloud in deep inelastic scattering. This is 
related to the pion contribution to the nucleon self-energy of Fig. [IJi. One needs to include 

d, but one also 
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the term in which the virtual photon interacts with the pion [18J, Fig. 
needs to include the effects of the virtual photon hitting the nucleon, Fig 
of momentum and charge would seem to require that the argument of the vertex function 
depends on the square of the invariant mass of the intermediate pion-baryon system (s) [T9] . 
Taking the form factor to have the standard form of depending on the square of the four- 
momentum transfer, between the initial nucleon and intermediate baryon (t), while natural, 
popular and effective [20j,[2j seemingly disagrees with charge and momentum conservation 
according to [T9] . 

But chiral symmetry (limit of vanishing pion mass) provides strong guidance. It is known 
that the ttN vertex function G n ^(t) and the nucleon axial form factor are related by the 
generalized Goldberger-Treiman relation |21j): 

MG A (t) = fA N (t), (!) 

where t is the square of the four-momentum transfered to the nucleons, Ga(£) is the axial 
vector form factor and /„- is the pion decay constant. The result Eq. ([!]), obtained from a 
matrix element of the axial vector current between two on-mass-shell nucleons, follows from 
PCAC and the pion pole dominance of the pseudoscalar current. Using Eq. ([!]) has obvious 
practical value because it relates an essentially unmeasurable quantity G n N with one Ga that 
is constrained by experiments. However the t dependence inherent in Eq. (JTJ) would seem to 
violate the purported consequence of momentum conservation. Similarly the pionic coupling 
between nucleons and A particles has an off-diagonal Goldberger-Treiman relation [22H21] . 
obtained using similar logic: 

2MC*(t) = fA NA (t), (2) 
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FIG. 1: (a) Pionic (dashed line) contribution to the nucleon (solid line) self-energy, (b) External 
interaction, x, with the pion (c) External interaction, x, with the intermediate nucleon. (d) Effect of 
27r-nucleon interaction. 

where is the Adler form factor [251 126], accessible in neutrino-nucleon interactions. 

The present manuscript develops a method that satisfies momentum conservation, utilizes 
Eq. ([I]) and involves only on-mass-shell nucleons. The key to removing ambiguities lies in 
evaluating the relevant Feynman diagrams by carrying out the integration over the four- 
momentum k by first integrating over k~ (the light front energy) in such a way that the 
intermediate baryon is projected onto its mass shell. This allows the use of the on-mass shell 
form factors Eqs. and is manifestly consistent with charge and momentum conservation. 

Consider the contribution to the nucleon self-energy E 7r (iV), involving an intermediate 
nucleon, Fig. ([lja.), given by Feynman rules as 




where M, /x are the nucleon and pion masses. The quantity P represents the nucleon mo- 
mentum and spin, (p, s), evaluated in the proton rest frame. We use the notation: G n N{t) = 
g nN F(t) = fG A (t), with G A (0) = 1.267 ±0.04, M = 0.939 GeV, / w = 92.6 MeV,g wN = 
G wN (0) = 13.2 with F(0) = I. The term F(k 2 ) represents the pion nucleon form factor. Its 
dependence on a single variable is justified only if the pionic vertex function appears between 




E^iV) = -i3g 2 nN u(P) 



d 4 k 7 5 (^- Jt + M) 7 5 

(2tt) 4 (k 2 -fi 2 + iejtjp - k) 2 — M 2 + ie) 



u(P)F 2 (k 2 ), 



(3) 



two on-mass-shell nucleons. In that case, one may use a dispersion relation: 



1 f°° 

F{k 2 ) = - dt' Im[F{t')}/(k 2 - 1'). (4) 

K J(3mn) 2 

Performing the spin average of Eq. ^ leads to the result 

T ( N) = 3 £n [ d A k F 2 (fc 2 ) kjp 

n{ ] M J i(2ny (k 2 - n 2 + ie)((p - kf - M 2 + ie)' [> 

We evaluate E 7r (A^) using light-front coordinates: ^ = k° ± k 3 , k 2 = k + k~ — k\. Thus 
E 7r (iV) = / dk+d 2 k ± J, with 

j = 1 1 / dk~F 2 (k 2 ) ^— ^ (61 

The expression Eq. (|4]) for F[k 2 ) is not written explicitly here because the analytic structure 
is the same as that of l/{k 2 — [i 2 + ie). If < k + < p + , the first pole in k~ is in the lower 
half k~ plane (LHP) (as are the ones arising from F(k 2 )) and the intermediate nucleon pole is 
in the upper half plane (UHP). We integrate over the UHP, so that the only pole we need to 
consider is the one in which the intermediate nucleon is on its mass shell and the momentum 
k is space-like. For k + < and k + > p + all of the poles are on the same side of the real axis, 
and one obtains 0. We take the residue of the integral for which the nucleon is on shell so 
that k~ = p~ — A p+^+ ■ Using the residue theorem and integrating over k + leads to the result 
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This result is obtained by using the pseuodscalar form of 7rN coupling in Eq. ([3]), but the use 
of pseudovector coupling would give the same result because the intermediate nucleon is on 
its mass shell. 

To proceed we use a specific form of the form factor F, the commonly used dipole 
parametrization 

F{Q 2 ) = 1/(1 + (Q 2 /M 2 A )) 2 , (8) 

with Ma as the so-called axial mass. The values of Ma are given by Ma = 1-03 ± 0.04 GeV 
as reviewed in [21]. This range is consistent with the one reported in a later review [27] . 
A somewhat lower value (0.85 Gev) is obtained [28j if one restricts the extraction region 
to very low values of Q 2 , but we need higher values to evaluate Eq. (J7|). Using this dipole 
parameterized form factor F gives 

Z n (N) = -3Mgl N -^— ^ x 

y7rAr 4(27T)3 6 (|_ 1 )V2 {a _ h y 
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[V(4 - 6)6 ( (a - 6) 2 (a(6 - 10) + 2(6 - 1)6) - 3a 2 (6 - 4) 2 61o£ 



+6 (4a 3 + a 2 (6 - 6)6((6 - 4)6 + 6) - 2a6 2 ((6 - 10)6 + 18) - 2(6 - 2)6 3 ) tan" 1 ^ - 1 
+6a6(6 - 4)V(a - 4)a(6 - 4)6tan- 1 ( J--1 )], a = ^ 2 /M 2 , b = M 2 A /M 2 . (9) 



To relate to chiral perturbation theory we expand in powers of a up to order /i 4 and 6 around 
unity to obtain a very accurate representation of the exact expression for < a < 0.04, 0.6 < 
6 < 1.6. We find 

EJN) = -3Mg 2 ^—[^^ + (-- - -^=) a + ixa^ 2 + 

2 104tt\ 2 16vra 8tt \ _ 2 Aog(a) 67tt 4V 

- + f= a ^ + f= I (6- 1) +a L ( 10 ) 

3 81^37 81^3 81^/ V 2 27^ 3^' 1 ; 

where the tilde indicates that a chiral expansion has been made. The term independent of 
the pion mass provides a -0.222 M correction to the bare nucleon mass, in contrast with 
an early approach (not using the heavy baryon expansion) which gives a contribution of 
formal order M(M/47r/ 7r ) 2 [29|. The term of order /i 3 reproduces the standard expression: 
-3g 2 A /(32nf 2 )^ M, ' 

The next step is to include terms with an intermediate A, the baryon excited state of lowest 
mass, which couples strongly to the ttN system. The effects of other intermediate baryons 
are not included in this first evaluation, but our technique can be applied to those states. We 
use the isospin-invariant interaction Lagrangian of the form C^na = S 2^A z (-p')g> J ' u u(p)d u 7r' 1 
+H.c. [221 [23] which yields the same result as the gauge invariant coupling of [12] for an on- 
shell intermediate A. We note that A* is a vector spinor in both spin and isospin space and 
g-wNA = v^6/2G 7r ATA(0), a notational relation between re-normalized coupling constants [22]. 
The contribution of the intermediate A to the nucleon self-energy is given by 

y (A) = i'?( 9 * NA Vr l (P) I — I + M A ) (g - kf 

* K 1 1 2M ] { ] J (2tt) 4 (k 2 - /i 2 + ze)((p - k) 2 - M 2 + ie) M\ 

xPfJ 2 \p-k)k^u(P)F 2 A (k 2 ), (11) 

where the factor of 2 arises from the isospin matrix element, Ma is the mass of the A and 
our notation for the projection operator P^J is given in [12] . We take the ratio of coupling 

is the S 
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constants to be { 9 g N ^ ) 2 = 72/25, which is the SU(Q) quark model result. The form factor F& 
is defined via G^na^) = g-n-N aFa^) = y^C^(t). Performing the spin average leads to the 
result 



S.(A) = 2(^) 2 i- ! d ^JMl { M 2 —p-k + MM A ) iP ~ k? 



2M M J z(2tt) 4 v f Ml 



2 r;2 {k-{ V -k))\ 1 



x-[k 



(p-k) 2 > (k 2 -/j 2 + ie)((p-k) 2 -M 2 A + i6) 



(12) 
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We evaluate S 7r (A) using light-front coordinates in a procedure analogous to that used for 
E„-(iV). The integral over k~ is done in the upper half k~ plane (UHP), so that the only pole 
is the one in which the intermediate A is on its mass shell and the momentum k is space-like. 
The result is 

QitNA-,2 11 1 f°° J4. F A(-t) ( \ , 1 i-mt-2. i\ /f2 , + \2 



^ V ; v 2M ; M(2vr) 3 3 7 (t + /j 2 ) \ 4Mf A J 

+ M A ) 2 + t) + Ml + ^{Ml — M 2 + t) 2 + 4tM 2 ) ■ (13) 

We turn to numerical evaluations. Lattice calculations [21] indicate that the ratio 
GnNA(t)/G n N(t) is constant as a function of the space-like values of t, thus here we use 



-^a(^) = F(t)- The integration of Eq. (13) yields a lengthy closed form expression. To gain 



insight, and compare with the general form of the chiral expansion of baryon masses in QCD 
e.g. [31H33] we take M A = M,b = l and expand in fx/M, (f - 0.72), f = Ml M ^ to find 

£.(A) = -™^?J^\\fM + (£ - 0-72)/ 2 (a)] (14) 

/ x (a) = -0.888a 2 + 1.01a 2 log(a) - 1.55a 2 (log(a) + 1.20) 

-0.402a 2 (log(a) + 1.24) - 0.00369a + 0.280a log(a) + 0.310 (15) 

f 2 (a) = (5.48a 2 + 1.46a 2 log(a) + 2.39a 2 (log(a) + 1.20) 

+0.128a 2 (log(a) + 1.24) + 1.02a + 0.318a log(a) - 0.0196), (16) 

where the tilde indicates that a chiral expansion has been made. The terms of order /z 4 log /i 2 
emphasized by [5],[3l] are included, but the expression also contains previously noted [33] 
dominating non-analytic terms of the form /i 2 log/i 2 . 

The total pionic contribution to the nucleon mass H n is given by 

£^EE^(iV) + ^(A), (17) 

and the chiral approximation is given by = E^iV) + E^A). These are shown in Fig. [i] 
as a function of the varying pion mass fi, the only parameter that is varied. Bare masses, 
M = 2.42 GeV, M = 2.06 GeV have been added to E^, T, n so as to reproduce the lattice 
data point at ^/47r/ w = 0.252(/x = 293) MeV. We use M A = 1.03 GeV. The use of the exact 
expression gives an approximately linear dependence on the pion mass, in agreement with the 
"surprisingly linear" results of lattice QCD simulations [H [35], found for values of \i greater 
than about 290 MeV. The LHP lattice data [35] are shown, and these are consistent with other 
lattice calculations as reviewed Varying the value of M A within the stated range changes the 
value of E„- only for \i > 0.5 GeV, and by 5 % or less. The low-order chiral approximation of 



Eq. (10) and Eq. (14) fails badly, showing that the chiral logarithms do not dominate for the 



relatively large values of \i used in many previous lattice QCD calculations. One could carry 



out the expansions of Eq. (10) and Eq. (14) to higher order in /i, but convergence requires 
many terms. One achieves a satisfactory description of S(iV) up to \x = 0.65 GeV by keeping 
terms up to order /i 24 , and of E(A) up to /i = Ma — M GeV by keeping terms up to order 



fi 20 . 
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FIG. 2: Nucleoli mass as a function of fi. Square blocks: LHP lattice data |35| . Solid: S w of Eq. (17), 
Eq. S and Eq. (Jl3|. Dashed: chiral approximation = E„.(JV) + E^A), Eq. (Fo| and Eq. (fuT 



It is worthwhile to compare our procedure with that of some others. For example, if one 
uses the heavy baryon limit to simplify Eq. (J3|, evaluates the integral by taking the pion to 
be on its mass shell and regularizes the divergent integral over momentum using a cutoff at 
a maximum momentum, one obtains results that correspond to the terms used in [5J. The 
relativistic procedure of |36j avoids the use of the heavy baryon limit by treating the nucleon 
recoil terms using an expansion procedure and uses dimensional regularization. We include 
all of the recoil terms and employ a cut-off procedure that is constrained by experimental 
data. In chiral perturbation theory, our procedure corresponds to keeping a specific set of 
higher-order terms with a fixed relation between them, a relation fixed by experimental data. 

Our results do not include contributions of order higher than l/f%. These may be con- 
sidered as keeping the lowest order pion cloud corrections using an expansion in powers of 
e = 1/(4tt^ j R) 2 , where R is a confinement radius PETJEB]. Here R ~ Vl2/M A , so e « 1/12. 
Thus we expect our results for the terms computed here to be accurate within about 10%. This 
argument was mainly applied to terms involving combinations of couplings of the nucleon to 
a single pion, but also holds for the n-pion-nucleon vertex e.g. as appearing in Fig. [l]i. These 
terms enter at higher orders in /i in chiral perturbation theory [39]. The coupling constant 
g n N and the confinement sizes of the pion and the nucleon, although not explicit in chiral 
perturbation theory, enter into the calculation of the diagram in terms of quarks and gluons 
and via the implicit dependence of f n and g n jy on the underlying strong coupling constant, 
as- Therefore we expect that the terms of the chiral Lagrangian will be consistent with the 
expansion in e. 

To test our treatment of the nucleon self-energy, we consider the contribution to lepton- 
nucleon deep inelastic scattering DIS arising from virtual pions. This is related to the term 
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Aijr, obtained from Feynman rules for the diagram of Fig lb, as 

8V, 

M* = 2M^. (18) 

This expression does not involve a "probability" , because the square of a nucleon light-front 
wave function does not appear. Note that charge and momentum are explicitly conserved: 
production of a pion of momentum k is accompanied by an intermediate nucleon of momentum 
p—k. The integrations over k~ , k± are carried out explicitly, and with the definition y = k + /p + 
one finds 

M„ = fdyfM, fM=f?{y)+£{v), 

fN( , _ , 2 i r t F\-t) 

U {y) { 2M ) (2vr)3 3 J {y , M , +y{Ml _ M , ))/{1 _ y) (t + /i 2 ) 2 

x (* + mi (iVf2 ~ Ml + 1)2 ) \ m + Ma)2 + 1] - (19) 



The functions f£ (y),f^(y) are shown in Fig. [3] where one observes that these functions are 
of roughly equal importance. 

The change in the quark distribution functions of the nucleon, 5qi(x), from this effect is 
given by the convolution formula as Sqi(x) = dyfTr(y)qJ(x/y), with qf the distribution 
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functions for quarks of flavor i in the pion. The related contribution to the nucleon structure 
function 5F 2 (x) is 



5F 2 (x) = f yfMFZ(x/y)dy, (20) 

J X 

where F^ is the pion structure function [HJ HD] • 

An integral involving the difference between the proton and neutron structure functions is 
particularly interesting: 

C 1 dr 2 r 1 

/ —(F*{z)-F?{z)) = 1/3-- dx(d(x)-u(x)), (21) 
Jo x ^ Jo 

where the first term, obtained if the bare nucleon has a symmetric sea, i.e. d = u, represents 
the Gottfried sum rule jH] . Experiment has clearly established violation of the Gottfried sum 
rule, and the most precise determination of the sea asymmetry [12] is 

D= [ (d(x) -u{x))dx = 0.118 ±0.012. (22) 
Jo 

Henley Sz Miller [10] showed that the pion cloud provides a natural explanation of the measured 
asymmetry. For Fig. the pion cloud of a proton will include n + (ud) and the 7r°, which 
has equal numbers of d and u. Only valence quarks of the pions are considered; the pion sea 
distributions are assumed to be symmetric. The probability for a 7r + n intermediate state is 
2/3, and that for a 7i°p state is 1/3. Including also the effects of an intermediate A leads to 

Dr = £dyv{lf?(v)-±tf{ y )), (23) 

with the probability of 7r~A ++ = 1/2 and that for 7r + A° =1/6. Since a bare baryon is assumed 
to have a symmetric sea, possible contributions from Fig. do not enter. Using = 1.03 
GeV, the nucleonic contribution is 0.173, and the A contribution is -0.064, so that the total 



is 0.109, within the experimental range of Eq. (22) 



To summarize: our light front treatment of the relevant Feynman diagrams reveals that 
the pion-baryon vertex function appears only between on-mass-shell baryons. This allows the 
vertex function to be expressed in terms of one variable, the invariant momentum transfer t, 
and to be constrained by experimental data. All ambiguities regarding the theoretical input 
needed to evaluate effects of the pion cloud to second-order in the coupling constants for 
the effects of intermediate N, A are resolved. The uncertainty due to the neglect of higher- 
order terms is estimated to be about 10%. Our procedure reproduces the observed linear 
dependence of the nucleon mass on the pion mass found in lattice QCD calculations and the 
flavor asymmetry of the nucleon sea. This work has implications for nucleon-nucleon scattering 
because one is instructed to use the coupling implied by Eq. (fTl), and also for computing pion 
cloud effects on the elastic electromagnetic form factors of nucleons [43J. 
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